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The  purpose  of  this  projeet  was  to  develop  the  modeling  and  analysis  of  both  PDE-based 
and  statistieal  point  proeess  models  for  repeat  aetivity  in  seeurity  applieations.  The 
researeh  team  has  made  great  strides  in  applying  sueh  ideas  to  urban  domestic  crime 
applications,  to  the  point  where  models  developed  in  research  papers  several  years  ago 
are  now  in  place  in  the  field  in  over  30  cities  worldwide.  This  project  developed  new 
basic  research  to  extend  many  of  these  ideas  beyond  domestic  crime  applications  to 
problems  abroad  involving  insurgents  and  also  to  other  areas  of  interest  to  the  Army  such 
as  email  traffic  on  social  networks  and  to  emotional  contagion  in  crowds.  The  research 
project  produced  15  journal  publications  and  an  additional  4  manuscripts  still  under 
review.  Published  and  submitted  journal  papers  included  the  work  of  two  postdoctoral 
scholars,  six  PhD  students  and  three  undergraduates,  in  addition  to  the  Pis.  For 
conciseness  we  report  on  those  manuscripts  that  were  published  or  submitted  in  the  last 
year  of  the  project.  Additional  manuscripts  and  publications  are  listed  for  this  award. 

Point  process  models: 

Crime  hotspot  maps  are  a  widely  used  and  successful  method  of  displaying  spatial  crime 
patterns  and  allocating  police  resources.  However,  hotspot  maps  are  often  created  over  a 
single  timescale  using  only  one  crime  type.  In  the  case  of  short-term  hotspot  maps  that 
utilize  several  weeks  of  crime  data,  risk  estimates  suffer  from  a  high  variance,  especially 
for  low  frequency  crimes  such  as  homicide.  Long-term  hotspot  maps  that  utilize  several 
years  of  data  fail  to  take  into  account  near-repeat  effects  and  emerging  hotspot  trends. 

We  show  how  point  process  models  of  crime  can  be  extended  to  include  leading  indicator 
crime  types,  while  capturing  both  short-term  and  long-term  patterns  of  risk,  through  a 
marked  point  process  approach.  Several  years  of  data  and  many  different  crime  types  are 
systematically  combined  to  yield  accurate  hotspot  maps  that  can  be  used  for  the  purpose 
of  predictive  policing  of  gun-related  crime.  We  apply  the  methodology  to  a  large,  open 
source  data  set  which  has  been  made  available  to  the  general  public  online  by  the 
Chicago  Police  Department.  This  work  was  published  in  the  International  Journal  of 
Forecasting. 

We  introduce  a  point  process  model  for  inter-gang  violence  driven  by  retaliation,  a  core 
feature  of  gang  behavior,  and  multi-party  inhibition.  Here,  a  coupled  system  of  state- 
dependent  jump  stochastic  differential  equations  is  used  to  model  the  conditional 
intensities  of  the  directed  network  of  gang  rivalries.  The  system  admits  an  exact 
simulation  strategy  based  upon  Poisson  thinning.  The  model  produces  a  wide  variety  of 
transient  or  stationary  weighted  network  configurations  and  we  investigate  under  what 
conditions  each  type  of  network  forms  in  the  continuum  limit.  We  then  fit  the  model  to 
gang  violence  data  provided  by  the  Hollenbeck  district  of  the  Los  Angeles  Police 
Department  to  measure  the  levels  of  excitation  and  inhibition  present  in  gang  violence 
dynamics,  as  well  as  the  stability  of  gang  rivalries  in  Hollenbeck.  This  work  was 
published  in  Discrete  and  Continuous  Dynamical  Systems. 

While  the  presence  of  clustering  in  crime  and  security  event  data  is  well  established,  the 
mechanisms  by  which  clustering  arises  is  not  fully  understood.  Both  contagion  models 
and  history  independent  correlation  models  are  applied,  but  not  simultaneously.  In  an 


attempt  to  disentagle  contagion  from  other  types  of  correlation,  we  consider  a  Hawkes 
process  with  background  rate  driven  by  a  log  Gaussian  Cox  process.  Our  inference 
methodology  is  an  efficient  Metropolis  adjusted  Langevin  algorithm  for  filtering  of  the 
intensity  and  estimation  of  the  model  parameters.  We  apply  the  methodology  to  property 
and  violent  crime  data  from  Chicago,  terrorist  attach  data  from  Northern  Ireland  and 
Israel,  and  civilian  casualty  data  from  Iraq.  For  each  data  set  we  quantify  the  uncertainty 
in  the  levels  of  contagion  vs.  history  independent  correlation.  This  work  was  published  in 
the  Annals  of  Applied  Statistics. 

Of  interest  to  ARO  is  the  analysis  of  patterns  in  e-mail  traffic.  Our  former  program 
manager  suggested  we  study  the  IKENET  dataset  from  West  Point.  Using  an  EM- 
type  approach,  we  fit  Hawkes  process  models  to  both  the  IKENET  and  the  Enron  e- 
mail  datasets.  We  show  that  the  self-exciting  models  adequately  capture  major 
temporal  clustering  features  in  the  data  and  perform  better  than  traditional 
stationary  Poisson  models.  We  also  investigate  how  accounting  for  diurnal  and 
weekly  trends  in  e-mail  activity  improves  the  overall  fit  to  the  observed  network 
data.  A  motivation  and  application  for  fitting  these  self-exciting  models  is  to  use 
parameter  estimates  to  characterize  important  e-mail  communication  behaviors 
such  as  the  baseline  sending  rates,  average  reply  rates,  and  average  response  times. 
A  primary  goal  is  to  use  these  features,  estimated  from  the  self-exciting  models,  to 
infer  the  underlying  leadership  status  of  users  in  the  West  Point  and  Enron 
networks.  This  work  is  under  review  at  J.  Amer.  Stat.  Assoc. 

Agent  based  and  PDE  models  (contagion,  Levy  flights,  and  optimal  control): 

We  consider  an  agent-based  model  of  emotional  contagion  coupled  with  motion  in  one 
dimension  that  has  recently  been  studied  in  the  computer  science  community.  The  model 
involves  movement  with  a  speed  proportional  to  a  “fear"  variable  that  undergoes  a 
temporal  consensus  averaging  based  on  distance  to  other  agents.  We  study  the  effect  of 
Riemann  initial  data  for  this  problem,  leading  to  shock  dynamics  that  are  studied  both 
within  the  agent-based  model  as  well  as  in  a  continuum  limit.  We  examine  the  behavior 
of  the  model  under  distinguished  limits  as  the  characteristic  contagion  interaction 
distance  and  the  interaction  timescale  both  approach  zero.  The  limiting  behavior  is 
related  to  a  classical  model  for  pressureless  gas  dynamics  with  “sticky"  particles.  In 
comparison,  we  observe  a  threshold  for  the  interaction  distance  vs.  interaction  timescale 
that  produce  qualitatively  different  behavior  for  the  system  -  in  one  case  particle  paths  do 
not  cross  and  there  is  a  natural  Eulerian  limit  involving  nonlocal  interactions  and  in  the 
other  case  particle  paths  can  cross  and  one  may  consider  only  a  kinetic  model  in  the 
continuum  limit.  This  work  has  been  accepted  for  publication  in  the  Journal  of  Statistical 
Physics. 

In  earlier  work  we  published  a  burglary  hotspot  model  [M.  B.  Short,  M.  R.  D’Orsogna, 

V.  B.  Pasour,  G.  E.  Tita,  P.  J.  Brantingham,  A.  E.  Bertozzi,  and  E.  B.  Chayes,  Math. 
Models  Methods  Appl.  Set,  18  (2008),  pp.  1249-1267]for  the  formation  of  hotspots  of 
criminal  activity.  In  new  work,  we  extend  the  earlier  model  to  incorporate  a  more  realistic 
model  of  human  locomotion.  The  movement  of  the  criminal  agents  follows  a  biased  Eevy 


flight  with  step  sizes  distributed  according  to  a  power-law  distribution.  The  biased 
Brownian  motion  of  the  original  model  is  then  derived  as  a  special  case.  Starting  with  an 
agent-based  model,  we  derive  its  continuum  limit.  This  consists  of  two  equations  and 
involves  the  fractional  Laplacian  operator.  A  numerical  method  based  on  the  fast  Fourier 
transform  is  used  to  simulate  the  continuum  model;  these  simulations  compare  favorably 
with  the  direct  numerical  simulations  of  the  agent-based  model.  A  Turing-type  analysis  is 
performed  to  estimate  how  the  instability  of  the  homogeneous  steady  state,  as  well  as  the 
expected  number  of  hotspots,  depends  on  the  system  parameters  and  especially  the 
exponent  of  the  underlying  power  law.  The  assumptions  of  the  underlying  agent-based 
model  naturally  lead  to  a  separation  of  scales  of  the  diffusion  coefficients  in  the 
continuum  limit.  Using  these  assumptions,  we  asymptotically  construct  the  leading-order 
profile  of  the  localized  hotspot  of  criminal  activity.  This  work  was  published  in  the  SIAM 
J.  of  Applied  Math  and  involved  three  undergraduate  researchers  as  part  of  the  project. 

In  the  M3  AS  paper  cited  above  (2008),  the  work  used  a  static  or  otherwise  suboptimal 
police  response  to  crime  hotspots.  In  new  work,  we  introduce  a  program  of  police 
response  to  hotspots  of  crime  in  which  the  police  adapt  dynamically  to  changing  crime 
patterns.  In  particular,  we  choose  their  deployment  to  solve  an  optimal  control  problem  at 
every  time.  This  gives  rise  to  a  free  boundary  problem  for  the  police  deployment's  spatial 
support.  We  present  an  efficient  algorithm  for  solving  this  problem  numerically  and  show 
that  police  presence  can  prompt  surprising  interactions  among  adjacent  hotspots.  This 
work  was  published  in  Discrete  and  Continuous  Dynamical  Systems. 

Density  Estimation  Problem 

This  work  brings  ideas  from  nonlocal  means  image  processing  to  bear  on  the  problem  of 
statistical  density  estimation  as  applied  to  spatial  mapping.  Given  a  discrete  sample  of 
event  locations,  we  wish  to  produce  a  probability  density  that  models  the  relative 
probability  of  events  occurring  in  a  spatial  domain.  Standard  density  estimation 
techniques  do  not  incorporate  priors  informed  by  spatial  data.  Such  methods  can  result  in 
assigning  significant  positive  probability  to  locations  where  events  cannot  realistically 
occur.  In  particular,  when  modeling  residential  burglaries,  standard  density  estimation 
can  predict  residential  burglaries  occurring  where  there  are  no  residences.  Incorporating 
the  spatial  data  can  inform  the  valid  region  for  the  density.  When  modeling  very  few 
events,  additional  priors  can  help  to  correctly  fill  in  the  gaps.  Learning  and  enforcing 
correlation  between  spatial  data  and  event  data  can  yield  better  estimates  from  fewer 
events.  We  propose  a  nonlocal  version  of  Maximum  Penalized  Likelihood  Estimation 
based  on  the  HI  Sobolev  seminorm  regularizer  that  computes  nonlocal  weights  from 
spatial  data  to  obtain  more  spatially  accurate  density  estimates. We  evaluate  this  method 
in  application  to  a  residential  burglary  data  set  from  San  Fernando  Valley  with  the 
nonlocal  weights  informed  by  housing  data  or  a  satellite  image.  This  work  involved  all 
three  Pis  and  UCLA  PhD  student  Joseph  Woodworth.  This  work  was  submitted  to 
Philosophical  Transactions  of  the  Royal  Society  A. 


